arXiv:1509.04896vl [astro-ph.HE] 16 Sep 2015 


Mon. Not. R. Astron. Soc. 000, ??-[6](0000) 


Printed 17 September 2015 (MN style file v2.2) 


The extension of variability properties in gamma-ray bursts to 
blazars 

Qingwen Wu\ Bing Zhang^’^’^, Wei-Hua Lei^, Yuan-Chuan Zou^, En-Wei Liang^’®, Xinwu Cao 

^School of Physics, Hashing University of Science and Technology, Whang 430074, China; Corresponding author. E-mail: qwwu@hust.edu.cn 
^Department of Physics and Astronomy, University of Nevada Las Vegas, NV 89154, USA 
^Department of Astronomy, School of Physics, Peking University, Beijing 100871, China 
‘^Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 

^Department of Physics and GXU-NAOC Center for Astrophysics and Space Sciences, Guangxi University, Nanning, 530004, China 
^Guangxi Key Laboratory for the Relativistic Astrophysics, Nanning 530004, China 

^ SHAO-XMU Joint Center for Astrophysics, Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China 
^Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, 210008 Nanjing, China 


ABSTRACT 

Both gamma-ray bursts (GRBs) and blazars have relativistic jets pointing at a small angle from 
our line of sight. Several recent studies suggested that these two kinds of sources may share 
similar jet physics. In this work, we explore the variability properties for GRBs and blazars 
as a whole. We find that the cotTelation between minimum variability timescale (MTS) and 
Lorentz factor, T, as found only in GRBs by Sonbas et al. can be extended to blazars with a 
joint correlation of MTS ex The same applies to the MTS ex correlation 

as found in GRBs, which can be well extended into blazars as well. These results provide 
further evidence that the jets in these two kinds of sources are similar despite of the very 
different mass scale of their central engines. Further investigations of the physical origin of 
these cotTelations are needed, which can shed light on the nature of the jet physics. 

Key words: gamma-ray burst: general - BL Lacertae objects: general - galaxies: jets - meth¬ 
ods: statistical 


1 INTRODUCTION 


Relativistic jets are common in black-hole (BH) systems, where 
the BH mass ranges from stellar mass (~ IOMq, Mq is solar 
mass, e.g., in BH X-ray binaries, XRBs, and gamma-ray bursts, 
GRBs) to supermassive scale e.g., in active galac¬ 

tic nuclei, AGNs). Blazars are an extreme type of AGNs in which 
the broad band emission is dominated by non-thermal radiation 
that produced in a relativistic jet pointing towards us. Based on 
the equivalent width (EW) of emission lines, blazars are divided 
into two sub-types, i.e., flat spectrum radio quasars (FSRQs) and 
BL Lacertae objects (BL Lacs): BL Lacs have weak or no emis¬ 
sion lines (EW ^ SA ), whereas FSRQs hav e strong emission 
lines (EW > 5A, e.g.. lUrrv & Padov^ll995h . where the differ¬ 
ent emission-line pr operties may be triggered by transition of ac- 
cretion modes (e.g.. iGhisellini & Celotti[r200ll : lwu & Caoll2008l : 
IXu et alj|2009l : IWu et al. 2013h . The spectral energy distribution 
(SED) of blazars generally include two broad humps in the log v- 
log(izE„) diagram: a low energy hump peaking at infrared to soft 
X-ray band that is produced by synchrotron radiation, and a high 
energy hump peaking at MeV to GeV band that is normally at¬ 
tributed to inverse Compton (IC) scattering. GRBs are the most en¬ 
ergetic explosions in the universe with a total isotropic-equivalent 


energy 


10 “ 


erg 


(e.g.. lMeszarosll 200 ^ : IZhanglllOlIL 


for reviews and the references therein), which may be caused by 
either collapses of massive stars or mergers of compact objects. 
The emission of a GRB is believed to originate from an ultra- 
relativistic jet with a Lorentz factor of a few hundreds. Such a 
high Lore ntz factor is evidenced from the “com pactness” argu¬ 
ment (e.g.. iFenimore et '^ ll99ilZwetal .11201 ih or the early af- 
tergl ow light curves that show the signal of fireball deceleratio n 
(e.g.. lSari & Piranlll999l : iLiang et ^120101 : iGhirlanda et al.ll2012h . 
The prompt emission of GRBs in the sub-MeV energy range nor¬ 
mally exhibits a non-thermal nature with a s moothly-joint bro ken 
power-law shape (i.e. the “Band function”. [Band et al.lll993h . A 
thermal component superposed on the non-thermal component was 
also observed in some bursts. The radiation mechanism of GRBs 
is still hotly debated (e.g., sy nchrotron, inverse Compton, thermal 
emission o r a mix of them, see lZhan g|2014l : lKumar & Zhang|2oT^ : 
|Peeill201^ . for recent reviews and references therein). 

It is believed that the jet physics (e.g., formation mechanism, 
energy dissipation, radiation process etc.) in BH systems should 
be intrinsically similar, which is supported by the s o-called “fun¬ 
damental plane” as foun d in XRBs and AG Ns (e.g.. lMerloni et all 
I2OO3I : Ipalcke et al.ll2004l : lDong et alj|2014h . Despite the large dis¬ 
crepancy of the masses of AGNs and GRBs, several universal cor¬ 
relations have been found among the two types of objects, strength- 
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ening this possibility. For example. IWang & Weil (1201 ih found that 
the GRB afterglows and BL Lacs may have the same radiation 
mechanism based on a similar relation in radio luminosity and 
sp ectral slope i n radio and optical band, which is further confirmed 
by IWang et ^ 1 2014 ) with two new similar relations of spectral 
properties. IWu et ahT 201 ih found a universal correlation between 
synchrotron luminosity and Doppler factor in GRBs and blazars, 
where they assumed the prompt emission of GR Bs i s dominated 
by syn chrotron emission. iNemmen et al.l ( 12012|) and IZhang et al.l 
1 I2OI3I) found the efficiency of energy dissipation of jets is similar 
in GRBs and blazars. A similar case w a s also mad e for X RBs and 
low-luminosity AGNs dMa et al.|l2014l) . ILvu et al.l ( l2014h tried to 
explore the possible unified radiation physics in GRBs and blazars. 

Both GRBs and blazars are highly variable sources. The typi¬ 
cal minim um variability timescale (MTS) for blazars is around one 
day (e.g.. IVovk & Neronovll201^ . where ultra-fast variability as 
short as 3-5 min at TeV energies is also found in some blazars (e.g., 

lAharonian et alj[2007l : I Albert et alfcoOTl i lAleksic et all2014h . The 

typical MTS of GRBs i s around 0.5 s with the shortest timescale on 
the order of 10 ms (e.g.. lMacLachlan et al.l2013l : lGolkhou & Butleil 
|2014 , where the MTS is much shorter than the overall duration 
(e.g., Tgo). The physical origin of the flux variation is unclear 
for both GRBs and blazars, where several models have been pro¬ 
posed to explain their temporal variability. The conventional mod- 
els such as the internal shock scenari o (e.g., KobaTOshi_etd .119971) 
and the photospheric scenario (e.g., iLazzati et al. 2OO9I) link the 
rapid variability of GRBs directly to the activity of the central 
engine. An alternative scenario envisages that the variability is 
triggered in the emitting region through local relativistic motion 
(e.g., local relativistic turbulence or reconnecti on events), which is 
not re l ated to the central e ngine directly (e.g. iNaravan & Kuma3 
l2009h . IZhang & YanI (1201 ll) suggested that the temporal variabil¬ 
ity of GRBs may include two variability components: a broad 
(slow) component related to the central engine activity and a narrow 
(fast) component associated with relativistic magne tic reconnection 
and/or turbulence (see e.g. IZhang & Zhandl2014r) . This scenario 
can also explain the ultra-fast variability (e. g., several minutes) as 
found in some blazars at TeV energies (e.g. . lAharonian et al.l2007l : 
[Albert et al.ll2007l:lAleksic et al.ll20l4. 


Sonbas et alj ( 12014b found that the MTS becomes shorter as 


the gamma-ray luminosity (or Lorentz factor) increases in GRBs. 
In this work, we try to extend the gamma-ray variability prop¬ 
erties in GRBs to blazars, which is aimed to explore the possi¬ 
ble universal variability properties in them and shed light on the 
physical mechanism that triggered the flux variations in relativistic 
jets. Throughout this work, we assume the following cosmology: 
Hq — 70 km s“^Mpc“^, Do = 0.3 and Da = 0.7. 


2 SAMPLE 


We select the GRB sample from ISonbas et al.l ( l2014b . where the 
correlations between the MTS, Lorentz factors and Gamma-ray 
lumino sities were explored. The MTS of GRBs in l^nbas et al.l 
( I2OI4I) was derived using a Harr wavelet technique, but the MTS 
was simply de fined in reference to a noise floor in the measure¬ 
me nt (see alsolMacLachlan et alJl2013l) . iGolkhou & Butlej ( |2014|) 
and iGolkhou et all ( 20151) analyzed the MTS from the structure- 
function method based on non-decimated Haar wavelets for a large 
sample of GRBs from Swift and Fermi. They reported MTS val¬ 
ues that are less dependent on the measureme n t noise level com- 
pared former works (e.g.. IWalker et al] I2OOOI : iMacLachlan et al.l 


I2OI3I : l^nbas et alj|2014l) . Since this method get rid of artificial 
effect due top the no ise floor, hereafter we ad opt the MTS mea¬ 
surem ents reported in iGolkhou & Butleil (l2014b andiGolkhou et al.l 


( I2OI5I) . The Lore ntz facto r s rep orted i 


Sonbas et al 


2014 


___ taken from Lit et al.l ( 1201^ and Ghir]andaet_alJ(l20]_2). 
iGhirlanda et al.l ( l2012l) made use of the lBlandford & McKeel h97^ 
self-similar deceleration solution, which led to typically a factor of 
2 times smaller L than those derived from the conventional method 
for a same source. Since at the peak time, the blast wave is not in 
the Blandford-McKee self-similar regime, we believe tha t the con¬ 
ventio nal method is more reasonable to estimate F (see lLii et al.l 
I2OI2L for more discussions). We therefore select the Lor entz fac¬ 
tors that are estimated fr om the conventional method (e.g.- lLii et al.l 
l2012l : lLiang et aljboisl) . We neglect the sources with only a upper 
or lower limits for the Lorentz factor, so that the sample is mostly 
composed of those GRBs whose F was measured usin g the conven¬ 
tional afterglow onset method (e.g. lSari & Piranll999l) . Four GRBs 
with Lorentz facto rs estimated from the GeV data or optical flashes 
are also selected dHascoet et al.ll201^ . where the GeV and opti¬ 
cal flashes are assumed to be emitted from the external shock. The 
isotropic mean gamma-ray luminosities of GRBs are selected from 
iLii et alJ ( l2012l) . which are derived from the total isotropic energy, 
Eiso, and burst duration, Tgo, with L.y = {1 -\- z)Eiso/Tgo. The 
redshift, MTS, Lorentz factor, and L.y for these GRBs are listed in 
Table 1. In total, we include 28 GRBs, with 22 having estimated 

Lorentz factor values. _ 

_ A sample of 45 blazars are selected from IVovk & NeronovI 

( l2013l) . where the MTS is derived using a structure-function ap¬ 
proach for the gamma-ray data from the Fermi telescope. The 
Doppler factor of blazars can be estimated from the timescale and 
amplitude of radio flares by assuming the variability timescale of 
the flares correspond to the light-travel time across the emission re¬ 
gion and that the intrinsic brightness t emperature of the so urce is 
limited to the equipartition value (e.g.. lHovatta et al1l2009l) . Com¬ 
bined with the apparent jet speed, it is possible to calculate the 
Lorentz factor and jet viewing angle for a given blazar. There are 20 
blazars with estimated Lorentz factor s using this method, which is 
selected from ISavolainen et al.l ( 1201OL for more details and the ref¬ 
erences therein). The total 0.1-100 GeV fc-corrected gamma-ray lu¬ 
minosity, L.y, is calculate d from the energy fl ux (Dy) and the fitted 
parameters as reported in iNolan et al.l J2OI2I) . The redshift, MTS, 
Lorentz factor and L-^ for blazars are reported in Table 2. 


3 RESULTS 

ISonbas et alj ( 12014l) found an anti-correlation between MTS and 
the Lorentz factor in GRBs, i.e. MTS oc ® for F > 225, 

and a shallow- plateau at smaller Lorentz fact ors. With the MTS val - 
ues derived bv lGolkhou & Butleil ( |2014|) and iGolkhou et alj ( l2015h . 
we find a similar MTS — F anti-correlation for the 22 GRBs with 
estimated Lorentz factors, with no evidence of the shallow-plateau 
(GRB 100621A is an outlier, see Figure 1). After excluding GRB 
100621A, a best straight-line fit to GRB data with errorfl in both 
F and MTS gives (solid line in Figure 1) 

MTS 

log --= -(4.8 ± 1.5) log F -f (11.1 ± 3.2), (1) 

1 ~\~ z. 

^ We adopt a linear fitting that considers the uncertainties in 
both X- and y- coordinates provided in NUMERICAL RECIPES 
IN FORTRAN in all the linear fittings in this work (fitxy.for, 
http://www.nrbook.eom/a/bookfpdf.php 1 . 
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Table 1. The GRB sample 


Name 

z 

log MTS 

r 

logL-, 

Refs. 

Name 

z 

log MTS 

r 

log Lj 

Refs 

060210 

3.91 

-0.59+0.07 

264tt 

51.97 

1,3 

090510 

0.90 

-2.58+0.09 

750^5° 

53.61 

2,5 

060607A 

3.08 

0.24+0.07 

296+f 

51.57 

1,3 

090618 

0.54 

-0.73+0.41 


52.31 

2 

061007 

1.26 

-1.15+0.08 

436t3 

52.50 

1,3 

090812 

2.45 

-0.58+0.07 

501 

52.27 

1,6 

061121 

1.31 

-1.28+0.07 

175tl 

51.88 

1,3 

090902B 

1.82 

-2.50+0.15 

750j;l50 

53.62 

2,5 

070318 

0.84 

0.43+0.07 

1431^ 

50.62 

1,3 

090926A 

2.11 

-2.03+0.06 

800t®®® 

53.87 

2,5 

071010B 

0.95 

-0.34+0.07 

209^4 

51.15 

1,3 

091003 

0.90 

-1.67+0.12 


51.93 

1 

071031 

2.69 

0.73+0.08 

133^^’' 

50.98 

1,3 

091029 

2.75 

-0.17+0.08 

221 

51.85 

1,6 

080319B 

0.94 

-1.68+0.11 


52.65 

1 

100414 

1.37 

-1.96+0.14 


52.75 

2 

080319C 

1.95 

-0.22+0.07 

228tg 

52.35 

1,3 

100621A 

0.54 

-0.09+0.07 

52 

51.04 

1,6 

080804 

2.21 

-0.66+0.23 


52.43 

2 

100728B 

2.11 

-1.25+0.08 

373 

51.89 

2,6 

080810 

3.35 

-1.64+0.18 

409+®^ 

52.08 

2,3 

100906A 

1.73 

-0.63+0.08 

369 

51.90 

1,6 

090102 

1.55 

-1.89+0.18 

440 

52.94 

2,4 

110205A 

2.22 

-0.70+0.08 

177 

51.85 

1,6 

090323 

3.59 

-1.65+0.23 


53.08 

2 

110213A 

1.46 

-0.71+0.04 

223 

51.52 

2,6 

090424 

0.54 

-1.11+0.07 

300j:^® 

51.08 

1,3 

130427A 

0.34 

-1.35+0.07 

■^^^-25 

53.43 

1,5 


References: 1) and 2) represent the references for MTS measurements that selected from lGolkfiou & Butled i2014l) an d lGolldiou et alj j2015h respec¬ 
tively. __ 

The Lorentz fa cto rs are selected from 3 ) ILu et ayj20l^l^')[Ghiilanda et alj i201 A tuse reported peak time of afterglow in this work and equation 1 of 
iLii et alJ2012h : Si lHascoet et alJ i201.'5h and 6t lLiang et alJ boish . 


Table 2. The blazar sample 


Name 

z 

log MTS 

r 

log+7 

Name 

z 

log MTS 

r 

log+7 

IJy 0138-097 

1.03 

6.03 


47.02 + 0.05 

3C 273 

0.16 

4.7 


46.31 + 0.01 

4C -H15.05 

0.41 

5.78 

9.9 

45.82 + 0.06 

3C 279 

0.54 

5.48 

20.9 

47.59 + 0.01 

PKS 0208-512 

1.00 

5.61 


47.54 + 0.02 

PKS 1437-1-398 

0.34 

6.25 


45.30 + 0.13 

3C 66A 

0.44 

5.10 


47.26 + 0.02 

PKS 1510-089 

0.36 

3.84 

20.7 

47.40 + 0.01 

AO 0235-1-164 

0.94 

5.94 

12.1 

47.95 + 0.02 

RGB 1534-H372 

0.14 

6.42 


44.46 + 0.09 

PKS 0336-019 

0.85 

6.02 

23.0 

46.96 + 0.04 

4C -1-38.41 

1.84 

4.81 

30.5 

48.72 + 0.00 

PKS 0420-01 

0.92 

5.37 

11.4 

47.61 + 0.02 

NRAO 530 

0.90 

4.90 

39.0 

47.36 + 0.02 

PKS 0440-00 

0.84 

4.95 


47.57 + 0.02 

4C-1-51.37 

1.38 

5.70 


47.82 + 0.01 

PKS 454-234 

1.00 

4.83 


48.13 + 0.01 

IJy 1749701 

0.77 

4.68 

7.5 

46.82 + 0.03 

4C-02.19 

2.29 

5.73 

16.2 

48.01 + 0.05 

OT081 

0.32 

5.47 


46.18 + 0.02 

PKS 0521-36 

0.06 

4.57 


44.68 + 0.02 

S5 1803-H784 

0.68 

4.95 

9.5 

47.07 + 0.01 

PKS 0537-286 

3.10 

6.21 


48.59 + 0.03 

B2 1811-1-31 

0.12 

6.26 


44.86 + 0.04 

PKS 0537-H441 

0.89 

6.04 


48.19 + 0.04 

4C -1-56.27 

0.66 

6.60 

37.8 

46.93 + 0.02 

PKS0716-H714 

0.30 

4.80 

10.3 

46.73 + 0.01 

PKS 1830-210 

2.51 

4.44 


49.35 + 0.00 

PKS 0735-H17 

0.42 

6.05 


46.58 + 0.02 

IJy 2005-489 

0.07 

6.48 


44.77 + 0.03 

PKS 0829-H046 

0.17 

6.06 


45.61 + 0.02 

PKS 2052-47 

1.49 

5.57 


48.28 + 0.01 

4C -H71.07 

2.17 

4.41 

28.0 

48.35 + 0.03 

4C -02.81 

1.29 

5.88 


47.17 + 0.05 

OJ 287 

0.31 

5.11 

15.4 

46.12 + 0.02 

S3 2141-H17 

0.21 

5.45 


46.02 + 0.01 

S4 0917-1-44 

2.19 

5.08 


48.67 + 0.01 

PKS 2155-304 

0.12 

6.69 


45.99 + 0.01 

4C -H55.17 

0.90 

5.50 


47.60 + 0.02 

BL Lac 

0.07 

4.96 

11.6 

45.10 + 0.01 

S4 0954-1-658 

0.37 

5.67 


45.98 + 0.03 

B3 2247-1-381 

0.12 

6.15 


44.68 + 0.07 

4C -H29.45 

0.72 

5.59 

25.1 

47.30 + 0.01 

3C 454.3 

0.86 

2.94 

19.9 

48.74 + 0.01 

4C -H21.35 

0.43 

3.64 

45.5 

47.44 + 0.01 







where a 50% error is adopted for P considering uncertainties of 
the model parameters in deriving them (e.g., one order of mag¬ 
nitude uncertainty in the density of ambient medium and radia¬ 
tive e fficiency of the blast wave, IWang et ^ 1201 Sl : Izhang et al.l 
l2015h and 0.2 dex is adopted for MTS considering that the MTS 
measurements are derived from different telescopes with different 
wavebands (e.g.,MTS derived from Swift at 15-350 keV maybe 
1-1 .5 times longer than that derived from Fermi at 8-1000 keV, 
see Rjolkhou et al.ll201a . for more details). We find that this anti¬ 
correlation as found in GRBs can be extended to blazars, where 
20 blazars roughly stay on the best-fit line of GRBs (solid line in 
Figure 1). Therefore, we further present the best straight-line fit for 
these 41 sources (20 blazars -l- 21 GRBs), 


log 


MTS 
1 + 2 


= -(4.7 ± 0.3) log F + (10.7 ± 0.7), 


( 2 ) 


where a typical uncertai nty of 0.4 dex and 0.2 dex is assigned 
for MTS of blazars (e.g.. lVovk & Neronovll2013h and GRBs (e.g., 
iGolkhou et alJ201.5h . respe ctively, and a 50% or 30% error is i ntro- 
duced for F f or GRBs (e.g., Lii et al.l2012LlSonbas et alj2014ll and 
blazars (e.g.. lHovatta et aij|2009l) . respectively. It is found that the 
global slope for blazars/GRBs is roughly consistent with that found 
in 2 1 GRBs only within u ncertainties. 

ISonbas et al.l ( l2014h found that MTS is also anti-correlated 
with gamma-ray luminosity, L^, in GRBs, where the relation also 
becomes flat when L~, < lO^^erg s ~^. With the MTS value s 
adopted from iGolkhou & Butlej ( |2014|) and iGolkhou et al.l ( 12015l) . 
we find a similar anti-correlation, where the best straight-line fit 
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Figure 1. The minimum variability timescale, MTS, vs Lorentz factor, T, 
for GRBs and blazars. The solid line represents the best linear fit for 21 
GRB (open circles, excluding a GRB 100621A that deviate the correlation 
evidently), which can roughly extended to blazars. The dotted line repre¬ 
sents the best fit for GRBs (excluding GRB 100621 A) and blazars in our 
sample. 
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Figure 2. MTS vs gamma-ray luminosity, Lj, for GRBs and blazars. The 
solid line represents the best linear fit for GRB sources, where blazars 
roughly stay on the best-fit line of GRBs. The dotted line represents the 
best fit for blazars (squares) and GRBs(circles). 


to GRB data with errors in both and MTS gives (solid line in 
Figure 2) 

MTS 

log--= -(1.0±0.1)logL^ + (51.2 ±3.8), (3) 

1 ± -2 

wher e a typical uncertain ty of 0.2 dex i s adopted for b oth MTS 
(e.g. jGol^ou et ^l2015l) and L-y (e.g.. iLii et alj|20l'3) for these 


GRBs. The hlazars also roughly stay on the extension of the best-fit 
line of GRBs (solid line in Figure 2). We also give the best straight- 
line fit to data of GRBs and blazars considering the uncertainties of 
both quantities, which gives 

MTS 

log —— = -(1.1 ± 0.1) log Lj -I- (56.4 ± 0.7), (4) 

L Z 

where the typical u ncertainty of 0.4 dex is adopted for the 
MTS of blazars fe.g.. lVovk & Neronovll201^ . Again the slope of 
GRBs/blazars is quite consistent with that derived in pure GRBs. 


4 CONCLUSION AND DISCUSSION 


GRBs and blazars are two populations of sources with the most 
powerful and relativistic jets in the universe. Compared with the 
accretion physics, the jet physics is less clear. The MTS provides 
information on the size and location of gamma-ray emitting re¬ 
gion, which may shed light on the physical process responsible 
for the variability (e.g., central engine activity, the Lorentz factor 
of the jet, and probably local relativistic reconnection/turbulence 
in the jet). In this work, we explore the possible universal vari¬ 
ability properties in GRBs and blazars by compiling 28 GRBs and 
45 blazars with measured MTS, gamma-ray luminosities and esti¬ 
mated Lorentz factors from the literature. We find the correlations 
of MTS — r and MTS — L-y as found in GRBs can be naturally 
extended to blazars, which not only provide further evidence on the 
similarities of the jets in GRBs and blazars but also give possible 
new clues to understand jet physics in these BH systems. 

The slope of MTS — T correlation is —4.7 ± 0.3 for GRBs 
and blazars, which is quite consistent with —4.8 ± 1.5 as found 
in pure GRB sample. Incidentally, the slope of MTS — L.y for 
the entire GRB-blazar sample (—1.1 ± 0.1) is also quite con¬ 
sistent with that found in GRBs only (—1.0 ± 0.1). The simi¬ 
lar slopes of the GRB/blazar sample and the pure GRB sample 
strengthen that the correlations as found in GRBs can be extended 
to blazars. These results provide new similarities on the jet physics 
in GRBs and blazars, where several other similar properties have 
been found in last several y ears(e.g., IWang & Weill201 l| :|Wu et al 
I2OI 1: iNemmen et al.ll201 j). It should be noted that Sonbas et al. 


(1201^ found a knee-like feature in both MTS — T and MTS — L. 
correlations, where the anti-correlations become flat for T < 225 or 
Lj < 10®^erg/s. This break is not evident in our correlations, with 
only GRB 100621A having a quite small Lorentz factor (T = 52) 
and evidently deviating from the MTS — T correlation (it how¬ 
ever follows the MTS—L.y correlation v ery well). This is because 
we adopted the MTS valu es provided bv iGolkhou & Butl^ l l2014h 
and iGolkhou et alj (l2015ll . who effectively removed the artificial 
MTS values near the threshold due to the nose floor. Their wavelet 
method also ensures a better homogeneity between the GRB and 
blazar samples, which strengthens our joint analysis. 

There is no evident MTS-T correlation in the subsample of 
blazars (correlation coefficient r = —0.15 with a probability 
of p =0.5). The MTS and L-y is also only mildly correlated 
(r = —0.51 and p = 5 x 10“®) in pure blazars. These corre¬ 
lations are not as tight as those found in pure GRBs (|r| > 0.8 
with p < 10““*), which prevent us from comparing the slopes 
of the blazar subsample with those of GRB subsample directly. 
The poor correlations among blazars may be intrinsic, but could 
be also caused by some uncertainties or biases in measuring T and 
MTS. For example, the F is derived by combining the measure¬ 
ments of the Doppler factor and apparent jet speed which make 
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use of observational data in different energy bands. The emission 
from different bands may originate from different emission regions, 
which may have different Lorentz factor, so that the estimated T 
may not be reliable. Also the typical MTS measured for blazars 
is of the order of several months, the inadequate temporal cover¬ 
age (e.g. with Fermi) may cause observational selection effects. 
Future data with better quality and more samples are needed to 
address whether blazar subsample indeed has (or does not have) 
significant correlations. In any case, the fact that the correlations 
derived from GRBs can be extended to blazars is intriguing, which 
may hint towards some common physics behind two types of jets. 


The variability of GRBs and blazars has been historically used 
to understand the well known “compactness” problem, and the min¬ 
imum bulk Lorentz factor s can be estimated if the MTS is used 
as one of the inputs (e.g., Dondi & Ghisellinil 1 19951 : IPiranI 1 19991 : 
iLithwick & Saril200lh . lSonbas et al.l ( l2014l) proposed that the anti¬ 
correlation of MTS-T as found in GRBs can be understood by the 
opacity effect, even if the slope is different from the 77 opacity 
argument (the gamma-ray luminosity is not included in their anal¬ 
ysis). The transparency condition for pair production is r-,.y < 1, 
which implies MTS oc where j3 is photon index of the 

spectrum. The high-ener gy index 0 range s from 2 to 3 with an av¬ 
erage value of 2.6 (e.g.,lTang et al.ll2015h in the Band function of 
GRBs (e.g.. lBand et alJl9^ 1 and ~ 2 — 3 for blazars as derived 
from the indices of electron distribution above the broken elec¬ 
tron Lorentz factor in modeling their mult i -wave band SEDs (e.g., 
P2 ~ 3 - 5. iKang et al.ll2014l) . IWu et all ( I2OI ih found a univer¬ 
sal correlation between synchrotron luminosity, Lsyn, and Doppler 


factor, T>, in GRBs and blazars (Lsyn oc D' 


, 3 . 1 ± 0.1 


), where the 


prompt emission of GRBs and the first hump of blazars are assumed 
to be dominated by synchrotron emission. Assuming gamma-ray 
luminosities in the second hump is normally more or less similar 
to Lsyn of the first hump in blazars and D 2 r for the rela¬ 
tivistic jet with viewing angle 6 —>■ 0 °, we find L-y oc p3 i±o i 
If this is the case, we find MTS oc for the typical value of 

P 2.6 in blazars and GRBs, which is roughly consistent with that 
found in GRBs and/or blazars (Eq. 2). It should be cautious that 
the opacity effect just give the lower limit for the Lorentz factor, 
which is energy dependent (e.g., depend on the photons at highest 
energy) and this gives the uncertainties in the above explanation. 
It is interesting to note that the derived Lorentz factors from the 
transparency condition in GRBs and blazars are normally just sev¬ 
eral t imes lower than tha t derived from other independent methods 
(e.g.. IZou & Piranlboid : l^ang et al.ll2012h . which suggests that 
above explanation based on transparency condition may be reason¬ 
able. 


There might be another interpretation to the relation. Based 
on the observational data, o ne has the Lorentz factor ratio 
EGRB/Tbiazar ~ 500/10 (e.g.. lHovatta et al.l2009l : lLu et alJ2012l) 
and black hole mass ratio McRB/Mbiazar IOMq/IO^Mq 
(e.g.. lWang et al.ir2004h between GRBs and blazars. Even though 
the physics behind this scaling relation is not fully explored, 
it suggests M oc i.e. a central engine with a smaller 

mass launches a faster jet. Since the observed MTS is MTS ex 
i?/(r^c) (r^c5/eng)/(r^c) ~ 5/eng, where ieng OC 2GM/c? 

is the minimum variability time scale of the c entral engine (e.g. 
iKobavashi et ^Il997l : lzhang & Meszarosll2004) . one has MTS oc 
M, and hence, MTS oc F”"' ’’ (Eq.(2)). 
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